Some anisotropic universes in the presence of imperfect 
fluid coupling with spatial curvature 



We consider Bianchi VI spacetime, which also can be reduced to Bianchi types VIo-V- 
III-I. We initially consider the most general form of the energy-momentum tensor which 
yields anisotropic stress and heat flow. We then derive an energy-momentum tensor that 
couples with the spatial curvature in a way so as to cancel out the terms that arise due to 
the spatial curvature in the evolution equations of the Einstein field equations. We obtain 
exact solutions for the universes indefinetly expanding with constant mean deceleration 
parameter. The solutions are beriefly discussed for each Bianchi type. The dynamics of the 
models and fluid are examined briefly, and the models that can approach to isotropy are 
determined. We conclude that even if the observed universe is almost isotropic, this does 
not necessarily imply the isotropy of the fluid (e.g., dark energy) affecting the evolution of 
the universe within the context of general relativity. 

Keywords Bianchi models • Imperfect fluid ■ Anisotropic dark energy • Isotropization ■ 
Heat flow • Spatial curvature 

1 Introduction 

The Bianchi universes form almost complete class of spatially homogenous but not necessar- 
ily isotropic relativistic cosmological models. They provide generalizations of the standard 
Fricdmann-Lcmaitre (FL) models, which are based on the spatially homogeneous and isotropic 
Robertson- Walker (RW) metrics [U [5] . Such models are of interest in cosmology in favor of 
constructing more realistic models than the FLRW models with maximally symmetric spatial 
geometry. Namely, although the observed universe seems to be almost isotropic on large scales, 
the early and/or very late universe could be anisotropic pQ. Additionally, the interest in such 
models was promoted in recent years due to the debate that going on the analysis and the in- 
terpretation of the Wilkinson Microwave Anisotropy Probe (WMAP) [3J HI [51 [B] data, whether 
they need a Bianchi type morphology to be explained successfully [71 [51 [HI EH (HI H2] ■ 

While the spatial isotropy of RW metrics exclude the anisotropic stress and/or energy flux, 
once the metric is generalized to Bianchi types, perfect fluid representantion of the energy- 
momentum tensor (EMT) can also be generalized in accordance with the considered Bianchi 
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type; e.g., the EMT can yield off-diagonal terms and anisotropic stress. However, the investi- 
gation of the anisotropic stresses has been widely neglected until recent years [13]. That might 
be because of the expectation that a dominance of an anisotropic stress would give rise to an 
anisotropic expansion whereas the observed universe appears to expand isotropically. How- 
ever, the anisotropic stresses have been considered by some authors. For instance, Barrow [14] 
analyzed the cosmological evolution of matter sources that possess anisotropic pressures (by as- 
suming the anisotropy is very small) such as electric and magnetic fields, collisionless relativistic 
particles, spatial curvature anisotropies, etc. Besides these known sources, since the release of 
WMAP data, the possible inflaton field and dark energy candidates that yield anisotropic pres- 
sure have been gaining increasing interest in favor of resolving the likely anomalies in WMAP 
data. Several cosmological models which introduce a source that yield anisotropic stress to 
break the spatial isotropy, either during inflation or in the late time acceleration of the universe, 
have been proposed, e.g., by the anisotropic inflaton field (say, vectorial fields) HSl HB] , by the 
anisotropic cosmological constant or dark energy HEl H3 HH1 UHl HSl I2H I2H 123 123 123 . However, 
even if it turns out that these anomalies are statisticaly insignificant and standard A cold dark 
matter model fits the WMAP data very well, we still can not rule out such energy sources. The 
reason being that anisotropic stresses do not necessarily promote anisotropy in the expansion, 
and there can even exist such anisotropic stresses that can cause universe to isotropize as Bianchi 
type I (or type V) universes in the presence of isotropic fluid ([22])- 

The spatial curvature is referred as a possible source of anisotropic stresses in general rela- 
tivistic cosmologies by Barrow [T3 120] • According to this idea, the three-curvature terms in the 
Einstein tensor can be moved to the other side of the Einstein equations and reinterpreted as an 
additional "effective" EMT, and then the "effective" EMT can behave as an anisotropic source. 
Thus, any Bianchi type, with anisotropic spatial curvature, universe containing an isotropic fluid, 
can behave like a Bianchi type I universe containing that fluid plus an additional anisotropic 
fluid (which obtained by the method mentioned above). On the other hand, one may also set 
an interesting class of solutions by following a similar idea. One may assume an EMT which 
compensates the three-curvature terms that moved to the other side of the Einstein equations 
and the total effective EMT becomes isotropic whether the metric possesses anisotropic spatial 
curvature or not. Hence, in that case, any Bianchi type, with anisotropic spatial curvature, 
universe containing this EMT, can behave like a Bianchi type I universe containing an isotropic 
fluid. Accordingly, one may assume a fluid that couples with the spatial curvature of the uni- 
verse in a way so as to vanish the terms arise due to the deviation from spatially flat geometry 
in the evolution equations. Once that is assumed any anisotropy in the spatial curvature will 
give rise to an anisotropic stress in the fluid. There are no a priori nor observational reasons 
to exclude such a fluid from being candidate for the inflaton field and/or dark energy, whose 
natures are not clear yet. Conversely, such an approach can provide interesting generalizations 
of the known candidates for inflaton field and/or dark energy. 

This study based on the above idea is carried out in the framework of Bianchi type VI 
spacetime, which can be reduced to the Bianchi types VIo , HI with anisotropic spatial curvature 
and types V, I with isotropic spatial curvature. Compared the Bianchi types I [221 12Z] and V 
[2U, Bianchi types III [23 [221 ffl EE] and VI [23 [321 [23 [33 1251 I3H] spacetimes have been 
studied rarely in literature of general relativistic cosmology. The basic reason being that these 
spacetimes do not cover RW spacetimes. However, they can still approach arbitrarily close to RW 
spacetimes. Thus, these spacetimes can still be studied for constructing realistic cosmological 
models. 
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We, initially, do not constraint the EMT and give full freedom to its components conveniently 
with the considered metric; hence, it yields anisotropic pressure and heat flow. We do not 
ignore the heat flow basically in favor of not compromising from generality. Heat flow can be 
important in the early stages of the universe and has been studied by several authors, e.g., 
[571 1551 1551 HOI IH]. Beside that, the heat flow can give rise to a temperature gradient along a 
certain axis of the universe. Hence, it would also be necessary to examine whether it tends to 
zero in the late universe in our models, otherwise the universe would not appear isotropic. 

We obtain exact solutions for the universes that exhibit volumetric expansion with constant 
mean deceleration parameter. Such solutions are not inconsistent with observations, can also be 
appraised approximately valid for slowly time varying deceleration parameter and have widely 
been considered (e.g., [231 EHl EZJ S21 S31 S31 HS1 SB])- Such cosmological models have also been 
studied in the presence of various anisotropic fluids by several authors recently. Sharif and 
Zubair investigated Bianchi VIo models in the presence of cosmological constant, anisotropic 
dark energy and electromagnetic field [23]. Yadav and Yadav investigated a Bianchi III model 
in the presence of anisotropic dark energy by assuming the expansion scalar is proportional to 
shear scalar [25] ■ Akarsu and Kilinc [3D] investigated LRS Bianchi I models in the presence 
of perfect fluid and minimally interacting anisotropic dark energy. The same authors |22) pre- 
sented Bianchi III models, which can isotropize as the Bianchi type I models in the presence 
of isotropic fluid, in the presence of single anisotropic fluid. They [3T] also discussed the de 
Sitter volumetric expansion in the presence of anisotropic fluid in detail within the framework 
of Bianchi I spacetime and presented two models in the presence of an anisotropic fluid that 
obtained by minimally altering the conventional vacuum energy. 

This paper is constructed as follows: In Sec. 2, the metric and some basic kincmatical 
equations are introduced. In Sec. 3, we discuss some basic dynamics of the models in the context 
of general relativity and derive an EMT that couples with the spatial curvature in a way so as to 
cancel out the terms that arise due to the spatial curvature in the evolution equations. In Sec. 
4, we obtain exact solutions for the universes that exhibit volumetric expansion with constant 
mean deceleration parameter. We briefly discuss the kinematics and dynamics of the models, 
and determine the models that can approach to isotropy. 

2 The metric and some basic kinematical equations 

For convenience, we use the natural units, i.e., 8nG = 1 and c = 1. The Ricci tensor is defined 
as the contraction of the Riemann tensor as follows: R^ u = R x ^ u \, where fj,, v and A run from 
to 3. Similarly, for the spatial section of the metric, the three-Ricci tensor is defined as follows: 
3 Rij = 3 R k ijk, where i, j and k run over spatial components from I to 3. 

We consider the spatially homogenous but anisotropic spacetime described by Bianchi type 
VI metric in the form 

ds 2 = -dt 2 + A{tfe- 2mz dx 2 + B(t) 2 e 2nz dy 2 + C(t) 2 dz 2 , (I) 

where m and n are real constants, and A, B and C are functions of the cosmic time t only. It can 
be observed that the metric (I) can be reduced to Bianchi type VIo (BVIo) when n = m ^ 0, 
Bianchi type III (Bill) when n ~ ^ m, Bianchi type V (BV) when n = —m ^ and Bianchi 
type I (BI) when n = m = 0. While BI covers the spatially flat RW spacetime, BV covers the 
spatially open RW spacetime. On the other hand, the underlying Lie algebra of the isometry 
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group of BVIo and Bill are completely different from those of BV and BI, and they do not 
cover RW spacetimes. However, they can evolve to a phase where they are observationally 
indistinguishable from RW spacetimes. 

Non-zero components of the three-Ricci tensors ( 3 Rj) and the three-Ricci scalar ( 3 i? = 3 Rl), 
which shall be used in the derivation of the energy-mometum tensor that would be considered 
in the solutions, of the metric (1) are as follows: 

3 D i m2 — mn 3 r>2 n 2 — mn 3 3 m 2 + n 2 

^1 = ~(J2 ' R 2 = (J2 aild ^3 = £^ ( 2 ) 



and 



3 „ 2m 2 - 2mn + 2n 2 



R= ^ • ( 3 ) 

It can be observed that the three-Ricci tensors are not all identical unless n = — m. Strictly 
speaking, the metric has an anisotropic spatial curvature when n ^ —m (Bill and BVIo cases), 
but has an isotropic spatial curvature when n = — m (BI and BV cases). The spatial curvature 
parameters, m and n, shall appear as the constants of the motion in the evolution equations 
(see Eqs. (14)-(16)). Accordingly, the flatness of the space cannot change as the universe 
expands independent of whether the universe undergoes an accelerated expansion period or not. 
However, an accelerated period of expansion (strictly speaking, in the z axis, i.e., C, for the 
metric we consider) can drive the three-Ricci tensors and three-Ricci scalar indistinguishably 
close to zero. In other words, the universe can evolve to a phase where it is observationally 
indistinguishable from spatially flat geometry. 

The directional Hubble parameters, which express the expansion rates of the universe in the 
directions of x, y and z respectively, can be defined as follows: 

H! = j, ff 2 = f and tf 3 = §, (4) 

where an overdot denotes d/dt. The mean Hubble parameter, which expresses the volumetric 
expansion rate of the universe, can also be given as 

i—l 

where V = ABC is the volume scale factor of the universe. 

Two other kinematical quantities of observational interest in cosmology are the anisotropy 
parameter of the expansion A and the dimensionless mean deceleration parameter q, which are 
defined as 



3 H 2 3 V H 



d ( 1 
~dt \H 



(7) 



where a 2 = \(Jij(j % i (aij is the shear tensor) is the shear scalar. A is the measure of the deviation 
from isotropic expansion and the universe expands isotropically if A = 0. The universe exhibits 
accelerating volumetric expansion if —1 < q < 0, decelerating volumetric expansion if q > 0, 
and exhibits constant-rate volumetric expansion if q = 0. 
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3 Some basic dynamics and derivation of the energy-momentum 
tensor 



We assume that the fluid four- velocity is comoving; thus it^ = i5q with u^u 11 = — 1, where 8^ 
is the unit four-tensor. The EMT of the fluid within the framework of the metric (1) can most 
generally be written in the following form; 
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where p is the energy density of the fluid; pi, pi and p 3 are the directional pressures on the x, 
y and z axes respectively, and h 3 is the energy flux density along the z axis. Because the fluid 
is comoving, energy flux due to the bulk motion of the fluid is not possible. However, heat flow 
can contribute to the energy flux, and then, h 3 may be interpreted as the heat flow along the z 
axis, see [471148] , 

Because the fluid is in the rest frame, anisotropy in its pressure can arise due to the anisotropy 
in its equation of state (EoS) parameter, since the energy density is a scalar quantity. Accord- 
ingly, one may assume that the directional pressures of the anisotropic fluid are each proportional 
to its energy density, i.e., 

(Pl,P2,P3) = (WlP,W 2 p,W 3 p), (9) 

where w\, w% and w 3 are the directional EoS parameters on the x, y and z axes respectively 
and they can be functions of the cosmic time t. However, for brevity, we handle with pressures 
rather than directional EoS parameters in the stream of solutions, though we shall deal with 
them while discussing the physical behaviors of the models. 
The EMT given in (8) may be parametrized as follows: 



( 



\ 



-P 




H 





P + Pi 








P+P2 





h 3 



P+P3 / 



(10) 



where p is the deviation- free part of the pressure (its meaning shall be specified below) , and p\ , 
Pi and P3 are the deviations from p on the x, y and z axes respectively, i.e., 



Pi = Pi - P, Pi = Pi - P and pz = p 3 - p. 
The Einstein field equations, in mixed components, can be written as 



_rpfj, 
- 1 v 1 



(11) 



(12) 



where is the Einstein tensor and R — is the Ricci scalar. 

In a comoving coordinate system, Einstein field equations (12) within the framework of the 
metric (1), in case of (10), lead to the following system of equations: 



AB AC BC m 2 -mn + n 2 
AB + AC + BC C 2 



(13) 
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(14) 
(15) 
(16) 
(17) 

The conservation of the EMT tensor given in (10), . v = 0, leads to the following two 
equations; 

in — 7x 

p + 3(p + p)H + piHi +p2H 2 +p3H 3 = h 3 , 

h + 3hH - m{p 3 - px) + n{p 3 - p 2 ) = 0. (18) 

In fact, from twice contracted Bianchi identity we have G^- v — for the Einstein tensor and 
thus considering (12) one may observe that Bianchi identity already assures the conservation 
law for the EMT given in (10) in general relativity. One may check that the solutions presented 
in the following sections satisfy both of these equations. 

The energy density of the fluid can be written in terms of the geometrical parameters, by 
using the constraint equation (13) and the definitions of A and H, as follows: 

P = 3^(l--j ^ (19) 

The sum of the directional pressures can also be written in terms of the geometrical param- 
eters, by using (13)-(16), (19) and the definitions of A, H and q, as follows: 

2 / _ 3 „ \ m 2 — mn + n 2 



3p + p 1 +p 2 +p 3 = 3H z ^-l + 2q--Aj + — . (20) 

The behavior of A, which appears in both (19) and (20), is crucial for deciding whether the 
models approach isotropy or not. According to the conditions given by Collins and Hawking 
[49] . the model approaches isotropy if as t — » oo 

h 3 

i) V -> oo, ii) A -> and iii) p > and ► 0. 

P 

The first two of these conditions are geometrical constrains and the spatial section of the metric 
approaches isotropy if both of them are satisfied. Considering the physical ingredient of the 
models, we shall say that the models approach isotropy when the last condition (iii) is satisfied 
as well as (i) and (ii). The condition ^ — > is taken because otherwise there would be a 
temperature gradient along the z axis and the universe would not appear isotropic. 

To express A explicitly we first obtain the differences between the directional Hubble pa- 
rameters by using the evolution equations (14)-(16); 



a y 



Xb+a + J (pb~Pa+ 3 R b b - 3 K) Vdt 



and a > b, (21) 
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where a and b run over spatial components and summation convention is not applied for a and b 
indices, and A3, A5 and A4 = A3 + A5 are real constants of integration. Now using this equation 
together with H and the three- Ricci tensors (2), in the definition of A we obtain the anisotropy 
parameter of the expansion explicitly as follows: 



A 



1 1 



9H 2 V 2 
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(i,b=l 



(pb 
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3 R b b - 3 R%) Vdt 



and a > b. 



(22) 



One can observe that the integral term in (22), contains both the differences of the deviations 
from p (i.e, pb — p a ) and the differences of the three-Ricci tensors (i.e., 3 R\ — 3 R%)- Thus, 
trivially, in the presence of an isotropic fluid (p± = pi = P3) and isotropic spatial curvature 
( 3 R\ = 3 i?2 = 3 Rl) the integral term vanishes and the anisotropy parameter of the expansion 
is reduced to the following well known simple form (e.g., see |50) ): 



A = 



K- 



H 2 V 2 



(23) 



where K 2 {= | J2 a 6=1 [^a+&] 2 anc ^ a > b) is constant. 

However, this is not the only option that permits the reduction of the general form of the 
anisotropy parameter (22) to this simple form (23). There is yet another option, which was first 
proposed by Akarsu and Kilinc [22) . but without unveiling its relation with spatial curvature. 
One may allow coupling of the spatial curvature with the fluid, i.e., the content of the integral 
term may be interpreted as pb + 3 Rb and p a + 3 R%- In convenience with this interpretation, we 
may assume that the deviations from p (i.e., p\ , p2 and ps) originate from the spatial curvature in 
a way so as to cancel out the terms arise due to the spatial curvature in the evolution equations 
(14)-(16), i.e., 
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and so that, using the three-Ricci tensors and three-scalar, we have 
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(25) 



Doing that, as would be seen in the following section, the terms arise due to the spatial curvature 
of the metric disappear in the evolution equations. Now, finally, substituting (25) in (10) we 
obtain an EMT in the following form: 
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(26) 



which also satisfies equations (18) that arise from the conservation law of the EMT. One can 
observe that if the spatial curvature of the universe is isotropic (i.e., when n — —to), then 
the EMT is also isotropic. This is the case for BI and BV. On the other hand, if the spatial 
curvature of the universe is anisotropic (i.e., when n ^ —to), then the EMT is also anisotropic. 
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This is the case for Bill and BVIo- One can also check that, in case of Bill (n = ^ m) and 
when heat flow is null (ft. 3 = 0), the EMT is reduced to a form which is equivalent to the one 
given by Akarsu and Kilinc 22 . 

Upon deriving the EMT, the meaning of the deviation-free part of the pressure (p) can also 
be specified. Substituting (25) in (20) we obtain 

p=(-l + 2q-^A)H 2 . (27) 

Note that, all terms that appear in this equation are dependent only on the volumetric evolution 
of the universe, due to the fact that A, q and H are parameters which are only dependable on 
the volume scale factor V and its time derivatives. Thus, p is the extracted part of the pi, pi 
and P3 , which is independent of the type of the metric as long as the volume of the universe is 
independent of the the spatial curvature parameters (m and n). 

Below we give exact models in the presence of the EMT (26) that derived according to (25), 
for the universe that expands with constant mean deceleration parameter. 



4 The models 

In a comoving coordinate system, the Einstein field equations (12) within the framework of the 
metric (1) and with the EMT (26), reduce to 

A B A C B C m 2 - ran + n 2 _ 

AB + AC + BC C~ 2 ~ P ' (28) 
B C BC 

B + C + BC=~ P > (29) 
A C AC 

A + C + AC = ~ P > (30) 

A B AB , N 

T+o + To=-P- (31) 



A B AB 
B 




- 77 = *3. (32) 



This system of the equations is not fully determined, because there are five linearly independent 
equations (28)-(32), but six variables (A, B, C, p, p and /13). Thus an extra constraint is 
necessary to close the system. Accordingly, we assume that the mean deceleration parameter 
q is constant. Such an assumption is not inconsistent with observations {271 142| and can also 
be appraised approximately valid for slowly time varying deceleration parameter |43j. One can 
exhaust all the possible constant q values for an indefinitely expanding universe when considering 
the exponential volumetric expansion (33) (q = —1) and the power-law volumetric expansion 
(34)(g>-l); 

V = cie 3fc * (33) 

and 

V = ai t 3 \ (34) 
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where c±, k and a\, I are positive constants. From (5) and (7), q = —1 for the exponential 
volumetric expansion and q — l^ 1 — 1 for the power-law volumetric expansion. Thus, the 
models with the exponential volumetric expansion and power-law volumetric expansion for I > 1 
exhibit accelerating volumetric expansion. The model for I = 1 exhibits constant-rate volumetric 
expansion, the models for I < 1 exhibit decelerating volumetric expansion. Thus, one may, 
phenomenologically, consider the anisotropic fluid we derived in the context of inflaton field or 
dark energy in the models with exponential expansion and power-law expansion for I > 1. 

It can be observed that, for the volumetric expansion laws given in (33) and (34), q and H 
are independent of the spatial curvature parameters (m and n), so do A and p from (23) and 
(27) respectively. Hence, the kinematical properties of the models arc independent of the spatial 
curvature parameters, i.e., A, B, C, H, A are independent of m and n. On the other hand, 
the properties of the fluid (p, pi, P2, Ps, h 3 ) and, trivially, spatial curvature properties of the 
models (R, 3 i?p are dependent on the curvature parameters. 

4.1 The models for exponential volumetric expansion 

Using the evolution equations (29-31) for the exponential volumetric expansion (33), we obtain 
the scale factors as follows: 

. C\ , f c4+<=5 r -3kt 

A = e kt 3 k e j (35) 

C2C3 

B = c 2 e kt +^- 3k \ (36) 
C = c 3 e kt +^ 3k \ (37) 

where C2 and C3 are positive constants and C4 and c 5 are real constants. 
The mean Hubble parameter is, 

H = k (38) 
and the directional Hubble parameters on the x, y and z axes are, respectively, 

H x =k + ( Ci + c 5 )e~ 3kt , (39) 

H 2 = k - c 5 e' 3kt , (40) 

H 3 = k-c 4 e- 3kt . (41) 

Using the directional and mean Hubble parameters in (6), the anisotropy of the expansion 
is obtained as follows: 

A= ^ + C :? +C5 V 6fct - (42 ) 

Using the scale factor on the z axis (37) in (2) and (3), the three-Ricci tensors and the 
three-Ricci scalar are obtained as follows: 

6 R\ = = e Akt 3 k e (43) 

g n , n 2 — inn .ou^a.-sn 

R 2 = 9 — e 3 " . ( 44 ) 

C3 
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C3 Z 



and 



p 3 = -3k 2 - (C 4 2 + C4 C 5 + C5 2 )e- 6fc * - ^ e -2 fe t-|^e— _ 



, _ 2m 2 - 2mn + 2n 2 oki 2 c 4 - 3k t 

6 R= - e 2kt 3— e . (46) 

C3 2 

Using the scale factors (35)-(37) in the constraint equation (28), the energy density of the 
fluid is obtained as follows: 

„,2 / 1 ?\ -rh m 2 —mn + n 2 _oi. t _2 3 p -3jst 

p = 3k 2 - (c 4 2 + c 4 c 5 + c 5 2 )e 6kt - e 2kt 3 ~ e . (47) 

C3 2 

The deviation- free pressure can be obtained by using (36) and (37) in (29); 

V = -3fc 2 - (c 4 2 + c 4 c 5 + c 5 2 )e- 6kt . (48) 

Using (37) and (48) in (26), the pressure of the fluid on the x, y and z axes are obtained as 
follows: 

Pl = -3k 2 ( C4 2 + C4 C5 + c 5 2 )e~ ekt + ^e- 2kt -i^- 3kt , (49) 

^3 

p 2 = -3k 2 (c 4 2 + c 4 c 5 + c 5 2 )e~ ekt + H^e' 2kt ~i^- 3kt , (50) 

2 £4 ^-3fct 

Finally, using the the scale factors in (32), we obtain the heat flow as follows: 

h 3 = n(c 4 - c 5 )e~ 3fct - m(2c 4 + c 5 )e~ 3fct . (52) 

We observe that all the parameters are finite at t = 0. While the volume of the universe 
increases continuously, the scale factors may behave non-monotonically for small t values, but 
eventually all expand monotonically and diverge as t — > 00. A is null throughout the history of 
the universe if C4 = c 5 = 0; otherwise, it decreases monotonically and tends to zero as t — > 00. 
Energy density of the fluid may behave non-monotonically for small t values, but eventually 
tends to 3fc 2 as t — > 00. The heat flow is null throughout the history of the universe in two 
cases; i) if C4 = c${m + n)/(n — 2m) and ii) if C5 = and m = n/2 simultaneously; otherwise, 
it decreases monotonically and tends to zero as t — > 00. We also need to examine the late time 
behavior of h 3 /p and one can observe that it also tends to zero as t — > 00. Above analysis shows 
that the model approaches isotropy. The three-Ricci tensors and three-Ricci scalar may behave 
non-monotonically for small t values, but eventually tend to zero as t — > 00. Accordingly, the 
universe eventually evolves to a phase where it is indistinguishably close to a flat and isotropic 
spatial geometry independent of the three-curvature parameters m and n. The behaviors of 
the directional pressures may differ from each other and may be non-monotonic for small t 
values, but eventually all of them tend to —3k 2 as t — > 00. Thus the directional EoS parameters 
(wi = Pi/ p), for large t values, are as follows: w\ = W2 = W3 = — 1. Hence, the fluid we 
considered isotropizes and exhibits a behavior indistinguishable from that of the conventional 
vacuum energy (p vac /Pva,c = — 1, which is mathematically equivalent to the cosmological constant 
A), in relatively later times of the universe for this model. 
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Above analysis can be carried out for particular BVIo, Bill, BV and BI cases by giving 
appropriate values to the spatial curvature parameters m and n. However, as mentioned above, 
the kinematical properties of the models are independent of the spatial curvature parameters 
(hence, the type of the metric), but the dynamical quantities of the fluid, the heat flow, and 
trivially the spatial curvature properties of the models do. Thus, below, we briefly mention only 
some properties of the models that are particular to the considered Bianchi type. 

The case for Bianchi type VIo (n = m ^ 0): In this model /13 is null if C4 = — 2c§. 3 R{ = 
3 i?2 are nuu J but 3 i?| and 3 R are identical and may behave non-monotonically for small t values, 
but tend to zero as t — > 00. The directional pressures on the x, y axes are identical, and are 
always higher than the one on the z axis; i.e., p\ — P2 > P3- They may behave non-monotonically 
for small t values, but eventually all of them tend to —3k 2 as t — > 00. 

The case for Bianchi type III (n = ^ m): One can check that if C5 = — 2c<4, then the 
heat flow is null and the model would be equivalent to the one given by Akarsu and Kilinc [22] 
for exponential volumetric expansion. 3 R 2 is null, but 3 R\ — 3 i?| = 3 i?/2 and they may behave 
non-monotonically for small t values, but eventually tend to zero as t — > 00. The directional 
pressures on the x, z axes are identical, and are always lower than the one on the y axis. They 
may behave non-monotonically for small t values, but eventually all of them tend to —3k 2 as 
t — >• 00. 

The case for Bianchi type V (n = m 7^ 0): This model is equivalent to the one obtained 
by Singh et al. [39] in the presence of a perfect fluid with heat flow for constant mean deceleration 
parameter with q = — 1. The heat flow is null if C4 = 0. All the three- Ricci tensors are identical 
and equal to one third of the three-Ricci scalar. Thus, the curvature of the space is isotropic. 
They may behave non-monotonically for small t values, but eventually tend to zero as t — > 00. 
The directional pressures on the x, y and z axes are identical, thus the fluid is isotropic. They 
may behave non-monotonically for small t values, but eventually all of them tend to —3k 2 as 
t — >• 00. 

The case for Bianchi type I (n = m = 0): This model is the simplest case and equivalent 
to the one given by Kumar and Singh [37] in the presence of a perfect fluid for constant mean 
deceleration parameter with q — — 1. The heat flow is null. The three-Ricci tensors and the 
three-Ricci scalar are all null. Energy density of the fluid increases monotonically and tends to 
3k 2 as t 00. The directional pressures on the x, y and z axes are all identical and increase 
monotonically as the universe evolves and tend to —3k 2 as t — » 00. One may observe that the 
fluid lies in the phantom region, i.e, p/p < —1, except in limit. 

4.2 The models for power-law volumetric expansion 

Using the evolution equations (29-31) for the power-law volumetric expansion (34), we obtain 
the scale factors as follows: 

Ol 7 °4+°5 .,1-31 , . 

A = _!-t< e — fcr* , (53) 
a 2 a 3 

B = a 2 t l e^ tl ~ 3 ' , (54) 
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C = a 3 t l e^T t , (55) 

where 02 and a 3 are positive constants and 0,4 and 0,5 are real constants. 
The mean Hubble parameter is, 

H=\ (56) 
and the directional Hubble parameters on the x, y and z axes are, respectively, 

H l = l - + (a A + a 5 )t- 31 , (57) 

H 2 = l --a 5 r 3L , (58) 

H 3 = l --a 4 t- 31 . (59) 

Using the directional and mean Hubble parameters in (6), the anisotropy of the expansion 
is obtained as follows: 

A = 2 a 4 2 + 0,4 q 5 + a 5 2 ^ 2 _ 6 , 
3 

Using the scale factor on the z axis (55) in (2) and (3), the three-Ricci tensors and the 
three-Ricci scalar are obtained as follows: 



o _i in — inn 91 2 °4 ji-31 

3 i?i = r 2l e~— * , (61) 

»3 



q „9 n 2 — mn 9 , 2a 4 ,1-3; 

3 i?2 = 5 — t- 2l e~— * , 62) 

a3 2 

3 = M 
3 a 3 2 

and 

o 2m 2 — 2mn + 2n 2 9 , .204.1-31 

3 i? = 9 r 2l e~~ * . (64) 

a 3 2 

Using the scale factors (53)-(55) in the constraint equation (28), the energy density of the 
fluid is obtained as follows: 

3l 2 , 9 fi , m 2 — mn + n 2 20^ i_ 3 ! 

P=-a~ K 2 + a 4 a 5 + a 5 )*~ 6 5 £- 2 ' e -— * . (65) 

The deviation- free pressure can be obtained by using (54) and (55) in (29); 

P = ~ 3 ^ 2 + 21 - (a 4 2 + a 4 c 5 + a 5 2 )t~ 61 . (66) 

Using (55) and (66) in (26), the pressure of the fluid on the x, y and z axes are obtained as 
follows: 

— 3^ 2 + 21 , 9 fi , 7J 2 9, 2a 4 1-31 

Pi = 75 (a 4 2 + a 4 C5 + a5 2 )t- 6 ' + — i- 2 'e^— 4 , (67) 

t 03 
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— 3/ 2 + 21 , 9 ,. a, ITT 2 o, 2 "4 f l-3i 

P2 = 75 a 4 2 + a4C5 + a 5 t~ 6 e^'-i* , (68) 

— 3Z 2 + 2Z . 9 ,, fi , I1W! 9, 2a 4 i_ 3 j 

P3 = 75 («4 2 + a 4 c 5 + a 5 )*~ 6 ' -t- 2i e -— * . (69) 

E «3 

Finally, using the scale factors in (32), we obtain the heat flow as follows: 

h 3 = n(a 4 - a 5 )t~ 31 - m(2a 4 + a 5 )r 31 (70) 

The initial time of the universe is t = 0, it starts with zero volume and expands indefinitely for 
all values of I. However, because the expansion anisotropy and the three- Ricci scalar contribute 
the energy density of the fluid negatively, we need to examine (65) to figure out which values of 
I are convenient for which times of the universe by applying the condition p > 0. Accordingly, if 
a 4 = d5 = 0; the models for I < 1 may represent the relatively earlier times of the universe, the 
models for I > 1 may represent the relatively later times of the universe and the model for I = 1 
may represent the entire history of the universe provided that 3ct3 2 > ra 2 — ran + n 2 . Similarly, 
if a 2 + a 4 a 5 + a 5 2 ^ 0; the models for < I < 1/3 may represent the relatively earlier times of 
the universe, 1/3 < Z < 1 may represent the intermediate times of the universe and I > 1 may 
represent the relatively later times of the universe. Finally the model for I — 1 may represent the 
relatively later times of the universe if 3a 3 2 > m 2 — ran + n 2 , otherwise the intermediate times. 
We are particularly interested in the models that could eventually evolve to a phase where they 
are indistinguishably close to a flat and isotropic spatial geometry. Hence, the models for I > 1 
and the model for I = 1 provided that 3a3 2 > to 2 — ran + n 2 arc of our interest, since these 
are the only models that can represent the late universe and can expand indefinitely without 
violating the positiveness condition on the energy density of the fluid. One can observe that for 
these models, expansion anisotropy is null throughout the history of the universe if a 4 = a$ — 
and otherwise, it is larger for small t values but decreases as t increases and tends to zero as 
t — > oo. Besides, h 3 /p — >• tends to zero as t — > oo. This analysis shows that the models for 
I > 1 and the model for I = 1 provided that 3a3 2 > to 2 — ran + n 2 approaches isotropy. All 
the three-Ricci tensors and the three- Ricci scalar behave parallel to each other and may behave 
non-monotonically for small t values but all eventually tend to zero as t — > oo. Accordingly, 
the universes that exhibit accelerating power- law volumetric expansion (I > 1) and the universe 
that exhibits constant-rate volumetric expansion (1 = 1, provided that 3a 3 2 > ra 2 — ran + n 2 ), 
eventually evolve to a phase where they are indistinguishably close to a flat and isotropic spatial 
geometry independent of the three-curvature parameters ra and n. 

For I > 1 , the energy density of the fluid may behave non-monotonically for relatively small 
t values, but p = 3l 2 /t 2 for large t values. The directional pressures may behave differently 
and non-monotonically for small t values, but px = p 2 = ps = (— 3/ 2 + 21) /t 2 for large t values. 
Thus the directional EoS parameters are as follows w\ = u>2 — W3 = — 1 + -Jy for large t values. 
Accordingly, in the late universe, independently of the spatial curvature parameters, the fluid 
isotropizes and evolves to an isotropic fluid with an EoS parameter that lies between —1 and 
— | according to the value of I > 1. 

For I = 1 under the condition 3a 3 2 > ra 2 — ran + n 2 , the energy density of the fluid may be- 
have non-monotonically for small t values, but p = 3 ° 3 ~™^2 t 2 in ~ n for large t values. Similarly, 
the directional pressures may behave differently and non-monotonically for small t values, but 
Pi = /° 3 ; p 2 — /a 3 anc | p 3 _ -i-mn/a.3 f Qr j ar g e £ values. Thus the directional 
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1 2 -I 2 

EoS parameters, for large t values, are as follows: w\ = — + 3 " a3 2 , w 2 = — + j^p- arid 
W3 ~ — 3 "™^ , where u = 1 — (to 2 — mn + n 2 )/3a 3 2 . Thus, the fluid does not isotropize 
even in the very late universe for 1 = 1. However, if a 3 3> to and «3 > n simultaneously, then 
w'i = W2 — W3 = — ^. Hence, also in this model, the fluid may appear to be isotropic. 

Note that, while the late time dynamics of the fluid do not exhibit dependence to the 
spatial curvature parameters in case of I > 1, they do in case of I — 1. Thus below, while 
consdering reduced types, we shall only mention I = 1 cases for the late time behaviors of the 
EoS parameters. 

The case for Bianchi type VIo (n = m ^ 0): In this model, heat flow is null if CI4 = — 205. 
3 R\ and 3 R 2 are null, but 3 i?| and the three-Ricci scalar are identical and may behave non- 
monotonically for small t values, but eventually tend to zero as t — > 00. The directional pressures 
on the x, y axes are identical, and are always higher than the one on the z axis; i.e., p\ = P2 > P3, 
both for Z > 1 and 1 = 1. 

The model for I = 1 evolves to a phase where it is indistinguishably close to a flat and 
isotropic spatial geometry provided that 3a^ 2 > to 2 . The energy density of the fluid may 

o 2 _ 2 

behave non-monotonically for small t values, but p = 3 3 2 t™ f° r relatively large t values. 
Similarly, the directional pressures may behave differently and non-monotonically for small t 

— \-\-m 2 la 2 ^ la 2 

values, but pi = p 2 = — " 2 /a3 and p 3 = "' 2 '— for large t values. Thus the directional 

1 2 -.2 

parameters, for large t values, are as follows: w\ = w 2 — — + 3^2 and W3 = — ^ — 3 ™ aa i , 
where u = 1 — m 2 /2>a 3 2 . Thus, the fluid does not isotropize even in the very late universe. 
However, if 0,3 3> to then the analysis for the general solution is valid once more. 

The case for Bianchi type III (n = ^ m): One can check that if 05 = — 2a4, then the 
heat flow is null and the model would be equivalent to the one given by Akarsu and Kilinc [22] 
for power-law volumetric expansion. 3 R 2 is null, but 3 R{ = 3 i?| = 3 i?/2 and they may behave 
non-monotonically for small t values, but eventually tend to zero as t — > 00. The directional 
pressures on the x, z axes are identical, and are always lower than the one on the y axis; 
Pi = P3 < P2 for both I > 1 and I = 1. 

The model for I = 1 evolves to a phase where it is indistinguishably close to a flat and 
isotropic spatial geometry provided that 3a3 2 > to 2 . The energy density of the fluid may 

o 2 _ 2 

behave non-monotonically for small t values, but p = ^ or l ar g e t values. Similarly, 

the directional pressures may behave differently and non-monotonically for small t values, but 

— l+n 2 la 2 — 1+m 2 la 2 

Pi = pz = — " 2 and p 2 = — " 2 /a3 for large t values. Thus the directional parameters, 

1 2 1 2 

for large t values, are as follows: i«i = w 3 = — h Q " 2 and w 2 = — t- + q m ■> , where 

u = 1 — TO 2 /3ci3 2 . Thus, the fluid does not isotropize even in the very late universe. However, 
if (13 ^> to then the preceding analysis is valid once more. 

The case for Bianchi type V (n = — m ^ 0): This model is equivalent to the one obtained 
by Singh et al. 39] in the presence of a perfect fluid with heat flow for constant mean deceleration 
parameter with q > — 1. The heat flow is null if 04 = 0. All the three-Ricci tensors are identical 
and equal to one third of the three-Ricci scalar. Thus, the curvature of the space is isotropic. 
They may behave non-monotonically for small t values, but eventually tend to zero as t — > 00. 
The directional pressures on the x, y and z axes are identical, thus the fluid is isotropic. 
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For I = 1 under the condition a^ 2 > to 2 , the energy density of the fluid may behave non- 
monotonically for small t values, but p = 3aa a 7 t 3 ™ for large t values. The pressure may behave 

— l+m 2 la 2 

non-monotonically for relatively small t values, but pi = P2 = P3 — — "a for large t values. 
Thus the EoS parameter, for large t values, is as follows: w± — w 2 — u>3 = — ^- + 3 ™ 2 , where 
u = 1 — m 2 /a3 2 . 

The case for Bianchi type I (n = m 7^ 0): This model is the simplest case and equivalent 
to the one given by Kumar and Singh |27j in the presence of a perfect fluid for constant mean 
deceleration parameter with q > — 1. The heat flow is null. The three- Ricci tensors and the 
three-Ricci scalar are all null. The models for I > 1/3 and the model for I = 1/3 under 
the condition 1/3 > a 2 + 0405 + a^ 2 satisfy all the conditions for approaching isotropy. For 
I > 1/3, the energy density of the fluid may behave non-monotonically for small t values, but 
p = 2l 2 /t 2 for large t values. The pressure may behave non-monotonically for small t values, but 
Pi = P2 = P3 — {—3l 2 + 21) /t 2 for large t values. Thus the EoS parameter, for large t values, is 
as follows: w\ = W2 = W3 = —1+37- According to this, the fluid evolves to a fluid with an EoS 
parameter that lies between —1 and — ^ according to the value of I > 1/3. For 1 = 1/3 under 
the condition 1/3 > a 2 + 0405 + a 5 2 , the energy density and the pressure of the fluid decreases 
monotonically. Thus the EoS parameter, for large t values, is as follows: W\ = W2 = if 3 — 1. 

5 Conclusion 

We consider the Bianchi type VI spacetime, which can be reduced to Bianchi types VLj, V, III 
and I. We derive an energy-momentum tensor under the assumption that the fluid is coupled 
with the spatial curvature in a way so as to cancel out the terms arise due to the spatial curvature 
in the evolution equations. This assumption results in an anisotropic pressure in the energy- 
momentum tensor when the spatial curvature is anisotropic and isotropic pressure when the 
spatial curvature is isotropic. The energy-momentum tensor we derived secures the evolution of 
the space, strictly speaking the kinematical properties of the universe, to be independent of the 
spatial curvature properties. We obtain the exact models for the universes that are expanding 
indefinitely with constant deceleration parameter. In the presence of the fluid we derived the 
following: all the models that exhibit accelerating volumetric expansion and the model that 
exhibits constant-rate volumetric expansion for a particular case, eventually evolve to a phase 
where they are indistinguishably close to a flat and isotropic spatial geometry independent of 
whether the fluid is isotropic or anisotropic and the universe possesses isotropic spatial curvature 
or not. According to this, even if the observed universe is almost isotropic, this does not 
necessarily imply the isotropy of the fluid affecting the expansion of the universe within the 
context of general relativity. In particular, we cannot rule out the possible anisotropic character 
of the inflaton field and dark energy, whose physical natures are not yet understood well. 
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